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9^ Abstract 

The Tevatron, where the top quark was discovered, and the currently functional Large Hadron 
Collider (LHC), with copiously produced top pairs, enable a detailed study of top-quark properties. 
In particular, they can be used to test the couplings of the top quark to gauge bosons. Several 
extensions of the standard model (SM) can give rise to anomalous couplings of the top quark to 
gauge bosons, in particular, the gluons. In this work we examine how top-quark polarization, 
which is predicted to be negligibly small in the SM, can be used to measure chromomagnetic 
and chromoelectric couplings of the top quark to gluons. We lay special emphasis on the use of 
angular distributions and asymmetries of charged leptons arising from top decay as measures of 
top polarization and hence of these anomalous couplings. Sensitivities that may be reached at the 
Tevatron and the LHC are obtained. 



a sudhansu.biswal@gmail.com 
b saurabh@prl.res.in 
c pankajs@kias.re.kr 



1 



I. INTRODUCTION 



The Tevatron has shut down its operations last year after 8.7 fb -1 of accumulated data. 
The first run of the Large Hadron Collider (LHC) with yfs = 7 TeV was already completed 
last year. In that run, LHC achieved 5.1 fb -1 of integrated luminosity. This year it has 
started at y/s = 8 TeV and has been projected to collect 15 fb -1 of data. After completing 
its run at \/s = 8 TeV, it is expected to start running at \/s = 14 TeV in 2014. With the 
standard model (SM) cross section for top-pair production at \/s = 14 TeV predicted to be 
around 830 pb, the LHC will provide ample opportunity to study top properties in detail. 

The top quark is the heaviest fundamental particle discovered so far with its mass m t = 
173.2 ± 0.9 GeV [1]. Mainly for this reason, it is considered to be a strong player in the 
determination of the mechanism of electroweak symmetry breaking (EWSB). The other 
consequence of its large mass is that its life time is very short and decays rapidly before any 
non-perturbative QCD effects can force it into a bound state. Thus, its spin information is 
preserved in terms of the differential distribution of its decay products. So by studying the 
kinematical distributions of top decay products, it is, in principle, possible to measure top 
polarization in any top production process. 

While already enough information about the top quark is available, which shows consis- 
tency with SM expectations, future runs at the LHC will enable more precise determination 
of its properties. The most recent experimental value of the top-pair production cross sec- 
tion at the Tevatron by CDF with 4.6 fb" 1 of data is a(ti) = 7.5 ± 0.31 (stat) ± 0.34 (syst) 
pb [2] for m t = 172.5 GeV and is consistent with the measurements from D0[3]. These mea- 
surements are in good agreement with the SM prediction of cr(tt)g^ LO = 7.08lol24ia27 pb 
for m t = 173 GeV [4]. The ti cross section has also been measured at the LHC, with a 
value of 161.9 ±2.5(stat)l5;o(syst) from CMS for an integrated luminosity of 2.3 fb _1 [5], and 
186 ± 13(stat)±20(syst)±7(lum) from ATLAS [6], for an integrated luminosity of 2.05 fb" 1 , 
in agreement with predictions of the SM. 

There seem to be hints of new physics from the study of top-pair production at the 
Tevatron in the forward-backward asymmetry of the top quark beyond the SM. Recent 
measurements by CDF [7] and D0[8] give a larger value for the asymmetry than predicted 
by the SM. 

Experiments at the Tevatron and the LHC have also produced results on top spin cor- 
relations [9-12]. LHC also has results on top polarization [13]. These are consistent with 
expectations from SM. Particularly, top polarization in the SM is predicted to be nearly 
vanishing at the LHC because the dominant contributions come from strong interactions, 
and are therefore parity conserving. Thus any deviation from zero would signal physics 
beyond SM. The errors are however still large, and new physics is not precluded. 

Top polarization and its usefulness in the study of new physics scenarios has been exten- 
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sively treated in the literature (for some recent papers in the context of hadron colliders, see 
[14-19]). For example, in Ref. [16], it was shown how top polarization could be utilized to 
probe the Z' couplings in the Little Higgs (LH) Model. In Ref. [17], the authors showed how 
top polarization may be used to determine the parameters of the Two Higgs Doublet Model 
(THDM) and Minimal Supersymmetric extension of Standard Model (MSSM). The effect 
of anomalous Wtb couplings on top polarization in single-top production has been studied 
in Ref. [18]. Probe of CP violation in single-top production using the polarization of top 
has been discussed in Ref. [19]. Refs. [20] suggest utilizing top polarization as a probe of 
models for the top forward-backward asymmetry observed at the Tevatron. 

In this work, we study top-pair production at the Tevatron and the LHC in the presence 
of anomalous gluon couplings to a ti pair. In particular, we examine the possibility of using 
top polarization and other kinematical observables constructed from top decay products 
in the laboratory frame to measure these anomalous couplings. Top chromomagnetic and 
chromoelectric couplings which we study here could arise in the SM or from new interactions 
at loop level. While the CP-conserving chromomagnetic coupling can arise in the SM at 
one-loop [21], the CP- violating chromoelectric coupling can only be generated at 3-loop 
level in the SM. Chromomagnetic and chromoelectric dipole moments of the top have been 
calculated at loop level in various new physics models such as MSSM [22], THDM [23], LH 
model [24] and in models with unparticles [25]. 

We calculate our observables at the Tevatron and at the LHC with centre-of-mass (cm) 
energies of 7 TeV (LHC7), 8 TeV (LHC8) and 14 TeV (LHC14). We also look at the 
sensitivities achieved in all these scenarios including statistical uncertainties with integrated 
luminosities of 8 fb" 1 at the Tevatron, 5 fb^ 1 at LHC7, 10 fb" 1 at LHC8 and 10 fb" 1 for the 
case of LHC14. 

Anomalous ttg couplings have been studied by several authors in the context of top-pair 
[26-32], top-pair plus jet [33] and single-top production [34] at hadron colliders. In Ref. [35], 
the author has used spin correlations in top-pair production at hadron colliders to probe 
chromomagnetic and chromoelectric dipole moments of top quarks. CP violation in top-pair 
production at hadron colliders including top chromoelectric couplings is studied in [36]. 

Apart from having a direct effect on top-pair production at hadron colliders, chromomag- 
netic and chromoelectric dipole couplings can have an indirect effect and modify the decay 
rate of b — > 57 at loop level [38, 39]. Using the measured branching ratio Br {b — > 57) [39], 
tight bounds on the chromomagnetic dipole coupling p were extracted, viz., 0.03 < p < 0.01. 

At the Tevatron and at the LHC, the dominant process of top production, viz., top-pair 
production, takes place through chirality-conserving QCD couplings in the SM. Thus, in 
the SM, the top polarization in top-pair production can only occur through the electroweak 
quark-antiquark annihilation into a virtual Z and is negligibly small. Any new physics in 
which new couplings to top are chiral can increase top polarization. The measurement of top 



3 



polarization is thus an important tool to study new physics in top-pair production. However, 
top polarization can only be measured through the distributions of its decay products. 
Hence, any new physics in top decay may contaminate the measurement of top polarization, 
and therefore of the new physics contribution in top production. Assuming only SM particles, 
any new physics in top decay can be parameterized in terms of anomalous tbW couplings as 



iVlLPL + fl R P R ) + (pt-Pb)u(f2LPL + f2 R P R ) 

mw 



(1) 



where in SM f XL = 1 and f\ R = f 2 L = f2 R = 0. Under the assumptions that (i) anomalous 
tbW couplings are small, (ii) the top is on-shell and (hi) t —> bW + is the only decay channel, 
it was shown in Refs. [40] that the charged-lepton angular distributions are independent of 
the anomalous tbW couplings. Thus, one can say that the charged-lepton angular distribu- 
tions are clean and uncontaminated probes of top polarization and thus of any new physics 
responsible for top production. 

For the above reasons, we choose, apart from top polarization, an asymmetry constructed 
out of the azimuthal distribution of charged leptons arising from top decay. 

In our work, we concentrate on the leptonic decay state arising from either t or i in 
top pair production. That is, we look at observables constructed from the charged lepton 
produced in t (t) decay, while the i (t), can decay into either a leptonic or a hadronic final 
state. Often we do not distinguish between observables related to t and those related to 
t. Thus measurements made for the top quark could also be made for the top antiquark 
in the process of top-pair production, and the results combined. However, in this case, 
information on CP violation would be lost. For the specific case of measurement of the 
CP-violating chromoelectric coupling, t and t observables have to be treated separately, and 
the corresponding partial cross sections appropriately added or subtracted. 

The rest of the paper is organized as follows. In the next section we discuss the formalism 
and the framework of our work. In Section III we discuss the application of the framework 
to the process of inclusive top-pair production at the Tevatron and the LHC, and present 
our results for the observables like top polarization, charged-lepton angular distributions 
and the lepton azimuthal asymmetry. Section IV deals with the statistical sensitivity of 
our observables to the anomalous couplings. The following section contains the conlcusions. 
The Appendix lists the production spin density matrix elements at the parton level for gg 
and qq initial states. 



II. THE FRAMEWORK 

We now describe the formalism underlying our analysis. 

We define the top-quark anomalous couplings to gluons including chromomagnetic and 
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chromoelectric dipole form factors by the ttg vertex 

m t 

where p and p' are the chromomagnetic and chromoelectric form factors respectively, q v 
is momentum of the gluon and m t is the mass of top quark. Of these, the p term is CP 
conserving, whereas the p' term is CP violating. We will treat the form factors p and p' as 
complex. In the SM, both p and p' are zero at tree level. 

For the calculation of the final-state charged-lepton distributions arising from either of 
t or i, we use the spin density matrix formalism. Since the top width of about 1.5 GeV is 
small compared to its mass, the narrow-width approximation (NWA) 
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p 2 — m 2 + im t T t 



Stf-mi). (3) 



can be utilized to factor the squared amplitude into production and decay parts as 

Pt ~ 



-— = 7lS(p 2 - m 2 ) y pXX , pAA' (4) 



A, A' 

where p xx ' and T xx ' are respectively the 2x2 top production and decay spin density matrices 
and A, A' = ± denote the sign of the top helicity. After phase space integration of p xx ' we 
get the resulting polarization density matrix a xx ' . The (1,1) and (2,2) diagonal elements 
of cr AA ' are the cross sections for the production of positive and negative helicity tops and 
cr tot = a ++ + o is the total cross section. 

Using Eq. (4) we can write the partial cross section in the parton cm frame as 



do 



E 



da xx ' f(T xx ') 



k , dcos9 t \pfPe 



32(2ir) 4 T t m t 

x dcos9 t d cos 9i d(f>£ EidEi dp'fy, (5) 



where the 6-quark energy integral is replaced by an integral over the invariant mass p 2 ^ of the 
W boson, its polar-angle integral is carried out using the Dirac delta function of Eq. (3), and 
the average over its azimuthal angle is denoted by the angular brackets. We obtain analytical 
expressions for the spin density matrix for top-pair production including the contributions 
of anomalous tig couplings to linear order at the parton level. These expressions for gg and 
qq initial states are listed separately in the Appendix. Use has been made of the analytic 
manipulation program FORM [41]. The expressions for the top-decay spin density matrix 
has been evaluated without linear approximation in anomalous tbW couplings in Ref. [18]. 
However, since we plan to work to linear order also in the tbW anomalous couplings, and 
evaluate observables dependent only on lepton angular variables, we need not include the 
dependence on tbW anomalous couplings. 
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III. THE TOP-PAIR PRODUCTION PROCESS 



We make use of the analytical expressions for the spin density matrix for ti production 
including anomalous ttg couplings to linear order listed in the Appendix. We find that 
at linear order, the real part of the coupling p and the imaginary part of the coupling p' 
give significant contributions to the diagonal elements of production density matrix, which 
are the ones which contribute to top polarization. The off-diagonal elements of the matrix 
get contributions from Imp and Imp', but not from the real parts of the anomalous cou- 
plings. The parton-level distributions are convoluted with parton distributions, which we 
do numerically. 

We neglect all fermion masses except that of the top and set Vf& = 1. For numerical 
calculations, we use the leading-order parton distribution function (PDF) set CTEQ6L [42] 
with a factorization scale of mt = 173.2 GeV. We also evaluate the strong coupling at the 
same scale, a s (m t ) = 0.1085. We make use of the following values for other parameters: 
M w = 80.403 GeV, the electromagnetic coupling a em {m z ) = 1/128 and sin 2 6 W = 0.23. 
We neglect the electroweak contributions in the production process. We take only one 
coupling to be non-zero at a time in the analysis except in Section IV where we show how 
simultaneous limits on two of anomalous ttg couplings may be obtained. In evaluating the 
angular distributions of charged lepton from top decay, we impose following acceptance cuts: 
p l T > 20 GeV and rapidity \rj\ < 2.5 on charged lepton. 



A. Top polarization 

The degree of longitudinal polarization P t of the top quark is given by 

with an analogous expression for the polarization of the top antiquark. 

In the SM, P t is predicted to be zero at tree level for top-pair production neglecting 
the contributions of s-channel 7, Z exchange in qq annihilation. We find that including 
non-vanishing anomalous ttg couplings, there can be non-zero top polarization asymmetry. 
In the expressions for the spin density matrix for top-pair production, we see that the 
contributions of Imp' have opposite signs in (1,1) and (2,2) elements and hence leads to 
non-zero Pt while the contributions of Rep have the same sign in these elements and thus 
do not contribute to top polarization. Thus top polarization can utilized to measure the 
coupling Imp' independently of all other anomalous ttg couplings. The diagonal elements 
of the density matrix for t also show that the polarization of t is the same as that of t, 
confirming that P t + Pj is indeed a measure of CP violation, proportional to the CP-odd 
coupling Imp'. P t is shown in Fig. 1 as a function of anomalous coupling Imp' in the linear 
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FIG. 1. The top polarization Pt in ti production at the Tevatron (bottom left), LHC7 (bottom 
right), LHC8 (top left) and LHC14 (top right) as a function of the anomalous ttg coupling Imp'. 
The grey band shows the 3a error interval in the SM without any pr cut. 



approximation, for the Tevatron and for the LHC7, LHC8 and LHC14. The grey bands in 
the figures denote the 3a statistical uncertainty in the measurement of Pt- The grey band 
is the thinnest for the 14 TeV LHC because of the largest cross section and therefore the 
smallest statistical error. 

We also study the effect of top-px cut on top polarization. In Fig. 1, we show top 
polarization for two different values of px cut i.e. pt > 100 GeV and px < 100 GeV. We 
find that for low-p^ tops, the top polarization is larger compared to high-p^ tops for the 
LHC while for the Tevatron, this observation is opposite. At Tevatron, high-p^ tops tend 
to have higher degree of polarization. 

We can understand the observation regarding the Tevatron as follows: At the Tevatron, 
the qq contribution dominates. In the diagonal elements of the spin density matrix for the 
gg-initiated contribution shown in the Appendix, the coefficient of Imp' is proportional to 
sin 2 9t = {VtIv 1 ) 2 - It is the Imp' which gives rise to the polarization, and so Pt is proportional 
to (p T ) 2 at the Tevatron, and it increases with transverse momentum. As for the LHC, the 
result is not so easy to see. 
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B. Angular distributions of the charged lepton 

Top polarization can be determined through the angular distribution of its decay prod- 
ucts. In the SM, the dominant decay mode is t — > bW + , with a branching ratio (BR) of 
0.998, with the W + subsequently decaying to £ + ue (semileptonic decay, BR 1/9 for each 
lepton) or ud, cs (hadronic decay, BR 2/3). The angular distribution of a decay product / 
for a top-quark ensemble has the form 



Here Of is the angle between fermion / and the top spin vector in the top rest frame and 
P t (defined in Eq. (6)) is the degree of polarization of the top-quark ensemble. 1^ is the 
partial decay width and Kf is the spin analyzing power of /. Obviously, a larger Kf makes / 
a more sensitive probe of the top spin. The charged lepton and the d quark are the best spin 
analyzers with = ng= 1, while K„ t = k u = —0.30 and Kb = —K\y+ = ~ 0.39, all k values 
being at tree level [43]. Thus the £ + or d have the largest probability of being emitted in 
the direction of the top spin and the least probability in the direction opposite to the spin. 
Since at the LHC, the lepton energy and momentum can be measured with high precision, 
we focus on leptonic decays of the top. 

To reconstruct the top-rest frame, one needs full information about top momentum. 
However, due to the missing neutrino, it is not possible to reconstruct completely and un- 
ambiguously the top longitudinal momentum and thus, this incomplete information may lead 
to large systematic errors. In this work, we focus on laboratory-frame angular distributions 
of the charged lepton and thus do not require full determination of the top momentum. In 
this sense, the observables we construct are more robust against systematic errors. Also, as 
mentioned earlier and shown in Refs. [40], the charged-lepton angular distribution in the lab 
frame is independent of any new physics in top decay and is thus a clean and uncontaminated 
probe of new physics in top production. 

We first obtain the angular distribution of the charged lepton in the parton cm frame, by 
integrating over the lepton energy, with limits given by rriyy < 2(p t -pi) < mf. This integral 
can be done analytically, giving the following expression for the differential cross section in 
the parton cm frame: 




(7) 



dcos9 t d cos 9i d(pg 32 T t m t (2%) 



da 1 1 




g 4 A xx ' \A( P 2 w )\ 2 dpl 



(8) 



where 



A* 




(l-r 2 ) 2 (l±cos^)(l T A)(l + 2r 2 ) , 



(9) 



24(l-/3 iC os^) 3 £ 2 




(10) 
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FIG. 2. The normalized polar-angle distribution of the charged lepton in tt production at the 
LHC7 (left), LHC8 (centre) and LHC14 (right) for the SM and with anomalous ttg couplings. 
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FIG. 3. The normalized polar-angle distribution of the charged lepton in tt production at the 
Tevatron for the SM and with anomalous ttg couplings. 



Here r = mw/mt and cos 9a is the angle between the top quark and the charged lepton in 
top decay in the parton cm frame, given by 

cos 9u = cos 9t cos 9e + sin 9 t sin 9t cos (fre, (11) 

where 9i and cj)g are the lepton polar and azimuthal angles. 

In the lab frame, we define the lepton polar angle w.r.t. either beam direction as the z 
axis and the azimuthal w.r.t. the top-production plane chosen as the x-z plane, with the 
convention that the x component of the top momentum is positive. At the LHC, which is a 
symmetric collider, it is not possible to define a positive sense for the z axis. Hence lepton 
angular distribution is symmetric under interchange of 9i and n — 9t as well as of 4>e and 

27T - (f) e . 

We first look at the polar-angle distribution of the charged lepton and the effect on it 
of anomalous ttg couplings. As can be seen from Fig. 2, where we plot the polar-angle 
distribution for LHC7, LHC8 and LHC14, the normalized distributions are insensitive to 
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anomalous ttg couplings. On the other hand, for the Tevatron, the polar-angle distribution 
are found to be somewhat sensitive as can be seen from Fig. 3. The sensitivity of polar-angle 
distributions on the anomalous ttg couplings have been studied in detail in Ref. [28] for the 
Tevatron, LHC7 and LHC14. Our results for these distributions agree with them. It is 
interesting to note that even though it is possible in principle to have a forward-backward 
asymmetric distribution at the Tevatron, the chromomagnetic and chromoelectric couplings 
in Eq. (2) do not generate an asymmetry. 




7C/2 71 371/2 7C/2 71 3n/2 7C/2 7t 37C/2 7C/2 jt 371/2 

<h 4>i <h 4>i 

FIG. 4. The normalized azimuthal distribution of the charged lepton in ti production at the 
Tevatron, LHC7, LHC8 and LHC14 for anomalous ttg coupling Rep (bottom) and Imp'. Also 
shown in each case is the distribution for the SM 



We next look at the contributions of anomalous couplings to the azimuthal distribution 
of the charged lepton. In Figs. 4 we show the normalized azimuthal distribution of the 
charged lepton in a linear approximation of the couplings for Tevatron, LHC7, LHC8 and 
LHC14 for different values of Rep and Imp' taken non-zero one at a time. We see that the 
curves for the couplings Rep and Imp' peak near <p£ = and <pe = 2n. 

In principle, it is possible to separate the dependence on the two couplings by taking the 
sum and difference of the distributions for t and t. The difference would be CP odd, and 
hence dependent only on Imp', whereas the sum would be CP even, depending only on Rep. 

We now discuss two angular asymmetries which would serve as a measure of the anoma- 
lous couplings. The first depends on the polar-angle distributions of the charged leptons 
from t and i, and the second one on the azimuthal distributions. 



10 



C. Charge Asymmetry 



We first look at a CP-violating asymmetry which is generated by the difference in the 
charged-lepton polar-angle distributions arising from the top and the antitop. We define the 
charge asymmetry 

. / /i \ 1 f c ° se ° , nfda + do" \ 

A ^ 9 °) = o 77TT / <*COS0 - , (12) 

2o-sm{vo) i-cose,, \dcos9 d cos 9 J 

where da ± /d cos 9 denote the differential cross sections for £ + and i~ production from t and 
i decay respectively, and ctsm(#o) is the cross section for either £ + or t~ production, with 
a cut-off of #o in the forward and backward directions of the lepton. It is obvious that for 
#o = 0, the numerator of Eq. (12) vanishes, because it measures the difference in the £ + and 
£~ production rates at all angles, which is zero from charge conservation. However, with a 
cut-off 6*0, A c h(6'o) can be non-zero, and is a measure of CP violation. It can be seen from 
the equations in the Appendix that A c h(0o) is proportional to Imp'. 

We plot in Fig. 5 the cross sections for charged leptons £ ± coming from decay of top/anti- 
top in top pair production as a function of cut-off angle 9 . We see from the Fig. 5 that 
the deviation in the cross section is relatively large in the range [tt/8, 37r/8]. To optimize 




It/8 It/4 3n/8 It/2 It/8 lt/4 3n/8 It/2 It/8 It/4 3lt/S lt/2 It/8 It/4 3lt/S It/2 

G 6 8q 6q 



FIG. 5. The cross section as a function of cut-off angle 9$ of the charged lepton in top-pair 
production at the Tevatron (bottom-let), LHC7 (bottom-right), LHC8 (top-left) and LHC14 (top- 
right) for anomalous ttg coupling Imp'. The SM cross section is also shown in each case. 



the charge asymmetry of lepton, we choose the cut-off angle #o to be ir/8 and evaluate the 
asymmetry as a function of Imp'. In Fig. 6, we plot the charge asymmetry of the lepton 
as defined in Eq. (12) as a function of Imp' for chosen value 7r/8 of 9 for Tevatron, LHC7, 
LHC8 and LHC14. 
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FIG. 6. The charge asymmetry A^Oq) as a function of charged lepton in top-pair production 
at the Tevatron (bottom-left), LHC7 (bottom-right), LHC8 (top-left) and LHC14 (top-right) for 
anomalous ttg coupling Imp' for 9q 
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We also study the effect of top-pr cuts on the lepton charge asymmetry A SM (9 ). From 
the top panel of Fig. 5, we see that at the LHC, keeping low px top/anti-top would help in 
enhancing Asm(#o) while at the Tevatron, the reverse is true. So, we put a cut on top/anti- 
top pt < 50 GeV at the LHC and px > 100 GeV at the Tevatron. We show the effects 
of these px cuts on charge asymmetry in top panel of Fig. 6. We find that though the 
statistical uncertainties increase due to the reduction in number of events, the asymmetry 
is increased enough times to compensate the reduction in events and thus results in the 
enhancement of the limits obtained by Asm(^o) on Imp'. 



D. Azimuthal Asymmetry 

As can be seen from Fig. 4, the curves are well separated at the peaks for the chosen 
values of the anomalous ttg couplings and are also well separated from the curve for the 
SM. We define an azimuthal asymmetry for the lepton to quantify these differences in the 
distributions by 

_ o-(cos0i > 0) - o-(cos0i < 0) 
* ~ a(cos 4>i > 0) + a(cos & < 0) ' [ ' 

where the denominator is the total cross section. This azimuthal asymmetry is in fact the 
"left-right asymmetry" of the charged lepton at the LHC defined with respect to the beam 
direction, with the right hemisphere defined as that in which the top momentum lies, and 
the left one being the opposite one. Plots of A$ as a function of the couplings are shown in 
Figs. 7 and 8 for Tevatron, LHC7, LHC8 and LHC14. 
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FIG. 7. The azimuthal asymmetry of the charged lepton in ti production at the Tevatron (1st 
column), LHC7 (2nd column), LHC8 (3rd column) and LHC14 (4th column) for different anomalous 
ttg couplings. 



From Fig. 4, we see that the azimuthal distribution of the decay charged lepton is more 
sensitive to Rep than to Imp'. Hence, we expect that the azimuthal asymmetry we construct 
in Eq. (13) would be a sensitive probe of Rep. This fact can indeed be seen from Figs. 7 
and 8 where the straight line for Rep is steeper than for Imp' implying a more significant 
contribution from the former. The reason we get straight lines for individual contributions to 
the asymmetry is that we are working in a linear approximation for the anomalous couplings. 

As mentioned earlier in the context of distributions, the dependence on the two couplings 
can be separated by choosing the sum and difference of the azimuthal asymmetries for t and 
t. The difference being CP odd, would be dependent only on Imp'. 

We also study the behavior of in the presence of cuts on the top transverse momentum. 
In the top panel of the Fig. 7, we show the behavior of as functions of Rep and Imp' 
with no cut on the top transverse momentum. In the middle and lower panel, we show A^ 
when we consider tops with < 100 GeV and 200 GeV respectively. Similarly in Fig. 8 
we consider high-p^ tops to evaluate the asymmetry. In the top, the middle and the lower 
panel of Fig. 8, we show A^ as functions of Rep and Imp' for top quarks with px > 100 
GeV, 200 GeV and 400 GeV respectively. 

We find that for high-p-^ tops, the azimuthal distribution is relatively more peaked than 
for low-pr ones. The reason for this is the (1 — f3 t cos 8 t e) 3 factor in the denominator of Eqs. 
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FIG. 8. The azimuthal asymmetry of the charged lepton in ti production at the Tevatron (bottom- 
left), LHC7 (bottom-right), LHC8 (top-left) and LHC14 (top-right) for different anomalous ttg 
couplings. 



(9) and (10). Thus, when (3 t is large, the distribution tends to peak near and 2tt. As a 
result asymmetry is larger for high-pr tops. This effect can be seen from Figs. 7 and 8 
where it can be easily seen that as the pt of the top is increased, the azimuthal asymmetry 
is larger. Hence we conclude that asymmetry constructed from high-pr tops would be more 
useful in constraining the anomalous top-gluon couplings. From the Fig. 8 we see that 
the coupling Rep is more sensitive at the Tevatron and LHC14. Though the value of the 
asymmetry also increases for lower-pr tops, the statistics is very low in that region and thus 
we do not gain in sensitivity. 



IV. SENSITIVITY ANALYSIS FOR ANOMALOUS ttg COUPLINGS 

We now study the statistical significance of the observables discussed in the previous 
sections to the anomalous ttg couplings at the Tevatron, LHC7, LHC8 and LHC14. For 
Tevatron, LHC7, LHC8 and LHC14, we assume integrated luminosities of 8 fb _1 , 5 fb _1 , 10 
fb , and 10 fb _1 respectively. To obtain the 3a limit on the anomalous ttg couplings from a 
measurement of an observable, we find those values of the couplings for which the observable 
deviates by 3a from its SM value. The statistical uncertainty at in the measurement of any 
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generic asymmetry Ai is given by 

where Af M is the asymmetry predicted in the SM and Af is the total number of events 
predicted in the SM. We apply this to the various asymmetries we have discussed. In 
case of the top polarization asymmetry, the limits are obtained on the assumption that the 
polarization can be measured with 100% accuracy. 





Pt 




A ch (6 = tt/8) 




Imp' 


Rep 


Imp' 


Imp' 


Tevatron 
LHC7 
LHC8 
LHC14 


[-9.75, 9.75] xlO" 3 
[-1.91, 1.91]xl0" 3 
[-1.13, 1.13]xl0" 3 
[-5.59, 5.59] xlO" 4 


[-8.90, 8.90] xlO" 3 
[-4.31, 4.31] xlO" 3 
[-2.09, 2.09] xlO" 3 
[-7.96, 7.96] xlO" 4 


[-7.88, 7.88] xlO" 2 
[-1.93, 1.93] xlO" 2 
[-1.05, 1.05]xl0" 2 
[-5.38, 5.38] xlO" 3 


[-3.98,3.98] x 10" 2 
[-6.42,6.42] x 10" 2 
[-4.56,4.56] x 10" 2 
[-3.53,3.53] x 10" 2 



TABLE I. Individual limits on anomalous couplings Rep and Imp' which may be obtained by the 
measurement of the observables at Tevatron, LHC7, LHC8 and LHC14 with integrated luminosities 
of 8 fb" 1 , 5 fb" 1 , 10 fb -1 , and 10 fb" 1 respectively. 



The 3<t limits on Rep and Imp' are given in Table I where we assume only one anomalous 
coupling to be non-zero at a time. In case of the lepton distributions, we take into account 
only one leptonic channel. Including other leptonic decays of the top would improve the 
limits further. 





Pt 


A4, 


A ch (6 = n/8) 




Imp' 


Rep 


Imp' 


Imp' 


Tevatron 
LHC7 
LHC8 
LHC14 


[-1.50, 1.50] xlO" 2 
[-1.37, 1.37] xlO" 3 
[-7.82, 7.82] xlO" 4 
[-3.51, 3.51] xlO" 4 


[-2.65, 2.65] xlO" 2 
[-8.51, 8.51] xlO" 3 
[-4.41, 4.41] xlO" 3 
[-1.85, 1.85] xlO" 3 


[-3.35, 3.35] xlO" 1 
[-2.26, 2.26] xlO" 2 
[-1.20, 1.20] xlO" 2 
[-6.09, 6.09] xlO" 3 


[-1.78,1.78] x 10" 2 
[-1.14,1.14] x 10" 2 
[-5.19,5.19] x 10" 3 



TABLE II. Individual limits on anomalous couplings Rep and Imp', with a cut pt <100 GeV 
on the top transverse momentum (for A^, we take pt <50 GeV), which may be obtained by the 
measurement of the observables at Tevatron, LHC7, LHC8 and LHC14 with integrated luminosities 
of 8 fb" 1 , 5 fb" 1 , 10 fb" 1 , and 10 fb" 1 respectively. 



In Table II, we give the 3a limits on Rep and Imp' applying a cut pr < 100 GeV on the 
top transverse momentum. From the table, we find that though the asymmetry increases 
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for Imp' with the cut, the limits on it do not change much because of the opposite effect of 
reduction in statistics. On the other hand, the limits on Rep actually worsen because the 
top-pr cut reduces the asymmetry for Rep. 





Pt 




A ch (6 = tt/8) 
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Imp' 


Tevatron 


[-6.79, 6.79] xlO" 3 
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[-3.13, 3.13]xl(T 3 
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[-2.19, 2.19]xl(T 3 


[-1.91, 1.91]xl0~ 3 
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[-9.45, 9.45] xl(T 4 


[-6.56, 6.56] xlO" 4 


[-1.00, 1.00] xlO" 2 





TABLE III. Individual limits on anomalous couplings Rep and Imp', with a cut px >100 GeV on 
the top transverse momentum, which may be obtained by the measurement of the observables at 
Tevatron, LHC7, LHC8 and LHC14 with integrated luminosities of 8 fb" 1 , 5 fb _1 , 10 fb~\ and 10 
fb" 1 respectively. 



In Table III, we give the 3a limits on Rep and Imp' applying a cut pr > 100 GeV. 
From the table, we find that with this cut, the limits are more stringent for Rep since the 
asymmetry A<j, for it is steeper as compared to the value without cuts. On the other hand, 
the limits on Imp' actually worsen because the top-p^ cut reduces the asymmetry for Rep. 

We also obtain simultaneous limits (taking both Rep and Imp' non-zero simultaneously) 
on these anomalous couplings that may be obtained by combining the measurements at 
Tevatron with LHC7, LHC8 and LHC14 . 



LHC14 -Tevatron 
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FIG. 9. The la (central region), 2a (middle region) and 3a (outer region) CL regions in the Rep- 
Imp' plane allowed by the combined measurement of two observables at a time. The left, centre 
and right plots correspond to measurements at the combinations Tevatron-LHC7, Tevatron-LHC8 
and Tevatron-LHC14 respectively. The x 2 values for la, 2a and 3a CL intervals are 2.30, 6.18 and 
11.83 respectively for 2 parameters in the fit. 
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For this, we perform a y 2 analysis to fit all the observables to within fa of statistical 
errors in the measurement of the observable. We define the following x 2 function 



where the sum runs over the n observables measured and / is the degree of the confidence 
interval. Pj's are the values of the observables obtained by taking both anomalous couplings 
non-zero (and is a function of the couplings Rep and Imp') and Oj's are the values of the 
observables obtained in the SM. a^s are the statistical fluctuations in the measurement of 
the observables, given in Eq. (14). 

In Fig. 9, we show the la, la and 3a regions in Rep-Imp' plane allowed by combined 
measurement of asymmetry at different experiments, taken two at a time. For this, 
in the \ 2 function of Eq. (15), we have combined the measurement at Tevatron with the 
measurements at LHC7, LHC8 and LHC14. From among the three combinations shown in 
Fig. 9, we find that the strongest simultaneous limits come from the combined measurements 
at Tevatron and LHC14, viz., ±0.006 on Rep and ±0.04 on Imp', at the 3a level. 

We now describe some other relevant work on the determination of the chromomagnetic 
and chromoelectric form factors of the top at hadron colliders. Some earlier work [26- 
28] made projections for possible limits on the couplings which would be obtained at the 
Tevatron and the LHC. In Ref. [28], the authors study polar-angle, transverse momentum 
and energy distributions of charged leptons coming from top decay at the Tevatron and LHC. 
Our results on polar-angle distributions of charged leptons are in agreement with theirs for 
nonzero Rep. They find 10-15% deviations from the SM in the angular distributions for 
values of p around 0.01 and of p' around 0.05. They also study lepton energy distributions, 
which bring in dependence on anomalous tbW couplings. Refs. [26, 27] proposed utilizing 
cross section measurements at the LHC and the Tevatron to put limits on the anomalous 
couplings. With the available data from Tevatron, Choudhury et al.[29] using a slightly 
different notation, put limits on the new physics scale A. They conclude that the cross 
section measurements at Tevatron would put a lower bound on A of about 7.4 TeV and 9 
TeV for p = ±1 respectively which in our notation would translate to p ~ [—1.94, 2.36] x 10 -2 
while at LHC7 the lower bound on A is 10 TeV which is equivalent to p < 1.75 x 10~ 2 . Hioki 
and Ohkuma [30] in their work, which they consider an update of [27], find that the Tevatron 
cross section results give bounds —.01 < p < .01 and 0.38 < p < 0.41 and |p'| < 0.12, of 
which, only the region around p = and p' = survive on using early LHC data. They 
also studied the effect on top pr, polar-angle distributions and invariant mass distributions 
from various combinations of p and p' in the range of 0.1 — 0.4 and found them to be 
give significant deviations from SM predictions. Hesari and Najafabadi [31] studied the 
fraction of the gg fusion contribution in tt production cross section at the Tevatron and 
at the LHC and concluded that this fraction is more sensitive at the Tevatron than at 




(15) 
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the LHC. From Tevatron results with integrated luminosity of 1 fb _1 , they quote limits of 
1.1 < p < 0.6,0.8 < p' < 0.8. For the full luminosity accumulated at the Tevatron, they 
project limits of 0.03 < p < 1.5 and 0.37 < p' < 0.37 while at the LHC7, the limits the 
expect are 0.04 < p < 0.98 and 0.15 < p' < 0.15. They also studied the charge asymmetry 
of the top at the LHC and found very loose bounds on p from it, and no sensitivity to p' . 

Other authors have considered possible limits on couplings from more detailed observa- 
tions. In Refs. [32], the authors considered probing the CP- violating chromoelectric dipole 
moment utilizing T-odd correlations constructed from jet and lepton momenta and found 
the expected limit to be \p'\ < 5 x 10~ 3 . Ref. [33] studies the chromoelectric coupling utiliz- 
ing various momentum correlations in ti and ti plus one-jet processes and derive the limit 
of |p'| > 0.35. In Ref. [34] Rizzo studied anomalous ttg couplings in single-top production 
at the Tevatron and at the LHC and concluded that the limits from this channel are about 
one order of magnitude smaller than those from the pair production processes. Effect of p 
and p' on spin correlations in ti has been studied in Ref. [35] and bounds are found to be 
-0.7 < p < 0.6 and -0.5 < p' < 0.5. 

In all above papers, authors have considered anomalous ttg couplings to be real. We 
consider both chromomagnetic and chromoelectric form factors to be complex. We found, 
however, that our observables get contributions only from the real part of p and the imagi- 
nary part of p'. The projections for the best limits on Rep from our observables are sometimes 
an order of magnitude better than those obtained from cross sections. In our analysis, we 
found top polarization and charge asymmetry of the charged lepton (both of which are CP 
odd) to be dependent on the imaginary part of the chromoelectric form factor. In Ref. [36], 
the author has considered chromomagnetic and chromoelectric form factors to be complex 
and construct various CP violating observables to obtain constraints on real and imaginary 
parts of p' . At the Tevatron with 30 fb _1 , the author obtained limits of 2.4 x 10~ 18 cm g s and 
1.1 x 10~ 18 cm g s on Rep' and Imp' respectively which in our units translate to 2.13 x 10~ 2 
and 9.77 x 10" 3 . At LHC14 with 150 fb~\ the limits obtained are 5.2 x 10" 20 cm g s and 
2.5 x 10 -20 cm g s on Rep' and Imp' respectively. In our units, limits are 4.62 x 10~ 3 and 
1.91 x 10~ 3 on Rep' and Imp' respectively. Considering that the luminosities we use for our 
limits are much lower, our limits are comparable to theirs. 

More recently the authors of Ref. [37] have obtained the upper bounds on chromomagnetic 
dipole moment (Rep) as 0.085 using the available data for the ti production at Tevatron 
and LHC for 7 TeV. They have predicted that the sensitivity to probe this coupling will 
be improved by a factor upto 4 by the boosted top measurements at 14 TeV LHC. On the 
other hand we predicted that measurement of total cross section, azimuthal asymmetry for 
leptons and top polarization asymmetry at the 7 TeV LHC can be sensitive to Rep and Imp' 
if these couplings were of the order of 10~ 2 — 10~ 3 . We also found similar improvements of 
sensitivity at 14 TeV LHC. 
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V. CONCLUSIONS 



We have investigated the sensitivity of the Tevatron, LHC7, LHC8 and LHC14 to the 
anomalous ttg couplings in top-pair production followed by semileptonic decay of the top. We 
derived analytical expressions for the spin density matrix for top-quark production including 
the contributions of both real and imaginary parts of anomalous ttg couplings. We work 
in the linear approximation of anomalous couplings. We find that only Rep and Imp' give 
significant contributions to the spin density matrix at linear order. It may be noted that 
Imp and Rep' do not appear in the observables we consider. This may be understood from 
the fact that the observables are even under naive time reversal T, viz., under change of 
sign of all the momenta and spins, without an interchange of initial and final states. In 
such a case, from the CPT theorem, the observables which are CP even can only arise from 
dispersive parts of form factors (in this case Rep) and the CP-odd observable can arise only 
from absorptive parts (in this case Imp'). 

Longitudinal top polarization can be utilized to separate the contribution of Imp' com- 
pletely independent of all other anomalous couplings whereas the total cross section can be 
used to separate the contribution of Rep. 

Since top polarization can be measured only through the differential distribution of its 
decay products, we also study the angular distributions of the charged lepton coming from 
the decay of the top. Charged-lepton momenta are measurable very accurately at the LHC 
and charged leptons have the best spin analyzing power. Also, charged-lepton angular 
distributions have been shown to be independent of any NP in top decay. We find that the 
polar-angle distribution is not very sensitive to the anomalous couplings. On the other hand, 
the normalized azimuthal distribution is found to be sensitive to the anomalous couplings. 
The azimuthal distribution peaks close to <p = and = 2ir, and the values at the peaks 
are quite sensitive to the magnitude and the sign of the anomalous couplings In order to 
quantify this difference and to be statistically more sensitive, we construct an inegrated 
azimuthal asymmetry from the azimuthal distribution of charged lepton. 

We study the effects of top transverse momentum cuts on top polarization and azimuthal 
distributions. We find that the top pt cut may enhance or reduce the top polarization 
depending whether we take a sample of low-pr or high-p-r tops. In our analysis, we observed 
that for the Tevatron an ensemble of high-p-r tops have higher degree of longitudinal top 
polarization as the function of imaginary part of anomalous chromoelectric coupling p'. 
Conversely, an ensemble of low-py tops reduce the top polarization at Tevatron for non-zero 
Imp'. On the other hand, in the case of the LHC, the observation is reversed. For high-p^ 
tops, top polarization is small and vice-versa. 

We consider two angular asymmetries, leptonic charge asymmetry A c h(#o) an d leptonic 
left-right asymmetry which serve as measures of the anomalous couplings. We find 
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that the A c ^(9q) is proportional to Imp'. The difference in the £ + and i~ cross sections is 
relatively large in the range 6 : [vr/8, 37r/8]. We choose the cut-off 6 to be it/8 to maximize 
the A c h(#o)- We furthermore study the effects of top pt cuts on the charge asymmetry and 
conclude that A c h(9o) is enhanced in the low top-p^ region at the LHC while the reverse is 
true for Tevatron. 

The effect of top pt cuts on the angular distribution of charged leptons can be easily 
understood through Eq. (8). For high-p^ tops, the angular distributions peak at extreme 
values leading to larger azimuthal asymmetries A^. We also study the effect of these cuts 
on the limits obtained by the measurement of top polarization and azimuthal asymmetry. 
We infer that the high-p^ tops give large azimuthal asymmetries and can thus give more 
stringent limits on the chromomagnetic top-gluon coupling as compared to low-p^ tops. 

In conclusion, we have shown that top polarization, and subsequent decay-lepton distri- 
butions can be used to obtain fairly stringent limits on chromagnetic and chromoelectric top 
couplings from the existing Tevatron data as well as data soon to be available after the 8 
TeV run of the LHC. The limits could be improved by the future runs of the LHC. 
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Appendix A: Spin Density Matrix for top/anti-top in top-pair production with 
anomalous ttg couplings 

In this Appendix, we present the spin density matrix elements a 1 - 7 for the top quark in 
the top pair production process. We include contributions of all anomalous couplings to 
linear order. 

Some general considerations can be used to anticipate the structure of the density ma- 
trix. Writing the density matrix as a sum of various contributions, ctsm from the SM, and 
CRep.Rep'.imp.imp', the respective contributions from Rep, Rep', Imp and Imp', 

° lJ = OSM + a Rcp + ^Rcp' + "hap + vLfJ ■ ( A1 ) 

Then, Hermiticity of the density matrix gives 

a ±± = a ±±* ; (A2) 
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implying that the diagonal matrix elements are real, and 



Re a ±T 
Im <r ±=F 



Re a T± 
— Im a T± . 



(A3) 



Thus, the only imaginary contributions come in the off-diagonal elements, changing sign 
under helicity flip. 

Let us now see what transformation under naive time reversal T tells us. Under T, o l i is 
transformed to cr 1 - 5 *. We note that p' terms are odd under CP, and therefore under T. Hence 
terms with Rep' would change sign under T. However, since T does not interchange initial 
and final states, the the contribution from Imp' does not change sign, as it arises from the 
absorptive part of some amplitude in the underlying theory. By the same argument, the 
Imp term changes sign under T, even though the corresponding interaction is T invariant. 
We thus have the relations: 



Thus, the diagonal density matrix elements can depend only on the couplings Rep and Imp'. 

The following relations arise from the parity transformation P, which flips the signs of 
the helicities, using the fact that the p' couplings are odd under P, and that there is an extra 
phase factor of —1 in the transformation of the off-diagonal elements: 



Now the above equations, together with the hermiticity relations in eq. A3 tell us that 



We conclude from the above that the diagonal matrix elements, which are all real, can 
only get contributions from the SM, and from Rep and Imp'. The off-diagonal elements 
have no contribution from the SM (since the SM amplitudes are real at tree level), and 




(A4) 




Im °^p> = Im a Rcp>- 



(A5) 



Re a SM = °> 

Re <rg = 0, 



(A6) 
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get a real contribution from Imp' and an imaginary contribution from Imp. Also, the Imp' 
contribution in the diagonal elements and the Imp contribution in the off-diagonal elements 
change sign under helicity flip. 

The density matrix a for the spin of the top anti-quark is obtained by changing the sign of 
the Imp and Imp' terms only in the off-diagonal element of the spin density matrix for the top 
quark. This can be seen from the fact that under the operation of CPT, where T is naive time 
reversal (reversing the sign of all spins and momenta, without interchange of initial and final 
states), the top spin density matrix elements would be transformed to the complex conjugates 
of the corresponding anti-top spin density matrix elements, with the helicity indices changing 
sign. However, this applies only to the real parts of couplings. Contributions containing 
imaginary parts of anomalous couplings, which arise from absorptive parts of amplitudes 
in an underlying theory, would change sign under this operation. Thus, because of CPT 
invariance, 

(AT) 
(A8) 

*w = -<l> = (A9) 



-±± 

°SM ~ 


~ "SM — 


_±± 

"SM ) 


-±± 
"Rep ~ 


■ "Rep — 


_±± 
"Rep; 


•±± — _ 
Imp' 


"imp' ~~ 


_±± 
"Imp' ) 




= — Im 


"imp' 


Re ^S' 


= -Re 


"imp': 



Im og, = Im = -Im (A10) 



Re ag, = -Re = -Re <r±* (All) 



AA' 



'Imp 
mp' 

We denote the spin density matrix for the top quark in the gg-initiated process as a 
and that for the gluon-gluon fusion process as cr^ . The labels (in subscript) s, t, u in a^' 
denote the s-, t- and u-channels contributions respectively whereas st, su and tu denote the 
interference between s- and t-channels, s- and u-channels and t- and u-channels respectively. 
The A and A' are top helicities and may take values ±. The total spin-density matrix 
elements a xx ' for top quark in tt-pair production in the parton cm frame (at parton level) 
are written as : 





= <r++ - 
qq 


r °gg,s + a gg,t 


~ 99,u 


+ cr ++ 
^ °gg,st 


+ a ++ 
1 gg,su 


^ °gg,tu 


(A12) 


a + - 


= (J tq r 






+ a + ~ 

1 99,u 


+ a gg,st 


+ <J + ~ 
1 gg,su 


+ °~gg,tu 


(A13) 


a- + 








+ a 99 + ,u 


+ a gg,st 


+ a~ + 
1 gg,su 


+ a gg,tu 


(A14) 


a 


= CJ qq ~ 


h a 99,s " 


V a 99,t 


a gg,u 


+ a gg,st 


a gg,su 


+ a gg,tu 


(A15) 



where, 



°tt = [ 16 Re P ± 8 lm P' & sin2 e * + ( 1 + $ cos2 Q t) + 4r t ] (A16) 

= — ^s 2 V§ sin 29 1 A M Imp ft + Imp'] (A17) 

<7±± = C s g 2 [8 Rep ± 8 Imp' (3 t cos 2 9 t + {l- f3 2 t cos 2 9 t )] (A18) 
2(7 
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Here 



°ff* = y^ 2 [-16Rep(l - A cos0 t ) =f 4Imp' {/3 t (8r t - 1) - cos^(l + 12r t ) 
+Pt cos 2 0* (1 + 16r t ) - 3f3 2 cos 3 4 + 2/3'f cos 4 t } + 1 + Ar t - 16r 2 
-/3 3 cos 4 - 8/3 2 r t cos 2 9 t + /3 3 cos 3 9 t - (3? cos 4 t ] (A20) 

of* = —s 2 VI [i Imp sin 9 t ft {± (1 + 8r t ) =f 2ft cos + 4r t ) 

± 3/3 2 cos 2 t T 2/3 3 cos 3 t } + Imp' sin 9 t {2 cos 9 t (3 t (1 + 6r t ) 
-3/3 2 cos 2 t + 2(3? cos 3 t - (1 + 12r t ) }} 

?} 9 % = [- 16Re / ? ( 1 + & cos 0*) T 4Imp' {/3 t (8r t - 1) + cos t (l + 12r t 
+/3 t cos 2 t (l + 16r t ) + 3/3 2 cos 3 0* + 2/3 3 cos 4 t } + 1 + 4r 4 - 16r t 2 
+/3 3 cos 9 t - 8/3 2 r t cos 2 t - /3 3 cos 3 t - /3 4 cos 4 t ] 

4Tn = — ^v 7 ! [i Imp sin 9 t t {=F (1 + 8r t ) =f 2ft cos t (l + 4r t ) 

=F 3/3 2 cos 2 t =F 2/3 3 cos 3 t } + Imp' sin 9 t {2 cos t /3 t (1 + 6r t ) 
+3/3 2 cos 2 t + 2/3 3 cos 3 t + (1 + I2r t ) }] 
^ st s 2 [4 Rep(-1 + 2fo cos 9 t ) ± 4 Imp cos 9 t {2 - A cos 

-1 _ R 2 ™<= 2 fl. _ R 3 ™<= fl. ^r, 2 fl. 1 



CT_t + = S 



99,st 

a 



gg,su 



99, su 
CT ±± 

gg,tu 



o ±T 

99, tu 



9s 2 ' 32 s 



24 (i-m 2 )' 



(A21) 



(A22) 



C st s 2 [4 Rep(-1 + 2ft cos 9 t ) ±4 Imp cos 6 t {2 - ftcos0 t ) 
1 - f3 2 cos 2 t - /3 3 cos t sin 2 t ] (A24) 

s 2 VI sin 9 t [±i Imp/3 2 (/3 t - cos 9 t ) - Imp'(l - f3 t cos t + 4r t )] (A25) 

. [4 Rep(l + 2/3 t cos0 t ) ±4 Imp'cos0t(2 + ft cos0 4 ) 

-1 + /3 2 cos 2 4 - /3 3 cos 9 t sin 2 t ] (A26) 
2C 

— - s 2 VI sin 9 t [±i Imp/3 2 (/3 f + cos 9 t ) - Imp'(l + f3 t cos t + 4r t )] (A27) 
m t 

C tu s 2 [±4 Imp' A {4r i (sin 2 t + cos 4 t ) - /3 2 (1 - 3 cos 2 9 t + cos 4 9 t )} 
+(3 2 {4r t sin 2 t + (3 2 t cos 2 4 (1 - /3 2 cos 2 9 t )} ] (A28) 
— - s 2 VI sin 29 1 fit [±i Imp /3 3 sin 2 t + Imp'{-/3, 2 sin 2 t + 2r t }l (A29) 



c« = t4; ^. = M; g «= ,. ,x2 ( A3 °) 



^4 3 ^4 

^ = 24 («-m 2 ) 2i ^ = 64 S (I - m 2 ) (A31) 

3 ^4 ^4 

^ = 64 I («-m 2 ) ; ^ = 192 (t-rn?) (fl - m 2 ) (A32) 

r t = ^; Pt = \/l-4^ (A33) 
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